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ABSTRACT 


The  spectral  domain  solution  for  the  wavelength  and 
characteristic  impedance  of  a  millimeter  wave  fin- line  was 
originally  published  by  Knorr  and  Shayda.  The  dispersion 
equations  were  subsequently  reformulated  by  Knorr  in  a  form 
more  suitable  for  numerical  computation. 

This  thesis  presents  a  reformulation  of  the  equations 
for  characteristic  impedance  for  the  same  purpose.  The 
equations  have  been  used  to  implement  a  computer  program, 
FINIMP .  The  program  runs  smoothly  without  the  overflow  and 
underflow  problems  experienced  by  Knorr  and  Shayda.  FINIMP 
data  is  compared  with  other  existing  data  and  good  agreement 
is  shown  to  establish  the  correctness  of  the  FINIMP  numer¬ 
ical  results. 
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I.  INTRODUCTION 

A.  BACKGROUND  AND  RELATED  WORK 

The  study  of  electromagnetic  energy  transmission  is  but 
one  important  area  in  microwave  and  millimeter-wave  engi¬ 
neering,  where  the  electromagnetic  waves  are  travelling 
through  some  transmission  medium,  which  provides  the  link 
between  the  transmitting  and  receiving  part  of  a  transmis¬ 
sion  system. 

In  recent  years,  fin- line  has  gained  in  importance  as  a 
transmission  medium  in  millimeter  wave  circuit  constructions 
[Ref.  1-4].  Fin-line  has  been  found  superior  to  micro¬ 
strip  at  millimeter  wavelengths  as  the  former  provides  eased 
production  tolerances,  low  dispersion,  broad  single  mode 
bandwidth,  moderate  attenuation,  better  compatibility  with 
hybrid  devices,  greater  freedom  from  radiation  and  higher 
mode  propagation,  combined  with  the  ability  to  construct 
simple  transitions  to  conventional  rectangular  waveguide. 

Figure  1.1  shows  a  3-dimensional  view  of  fin- line.  The 
stucture  may  be  viewed  as  a  slotline  with  a  shield,  a  ridged 
waveguide  with  dielectric  or  a  slab  loaded  waveguide  with 
fins  . 

The  fin-line  structure  was  first  proposed  for  millimeter 
wave  integrated  circuits  in  1974  by  Meier  [Ref.  1].  An 
early  paper  by  Meier  described  the  propagation  mode  as  a 
variation  of  dominant  mode  in  ridged  waveguide.  His  proce¬ 
dure  requires  a  test  measurement  to  determine  the  equivalent 
dielectric  constant  of  the  fin- line  structure.  This  is  both 
expensive  and  time  consuming.  Knorr  and  Kuchler  [Ref.  5] 
presented  a  frequency  dependent  hybrid-mode  analysis  of  slot 
line  with  open  boundary  using  the  spectral  domain  technique 
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which  was  suggested  by  Itoh  and  Mittra  [Ref.  6],  in  1975. 
Subsequently,  the  paper  of  prime  importance  in  the  estab¬ 
lishment  of  spectral  domain  technique  for  analyzing  the 
fin-line  structure,  so-called  shielded  slotline,  was 
pesented  by  Knorr  and  Shayda  in  1980  [Ref.  3].  One  of  the 
advantages  of  the  spectral  domain  approach  is  that  it  is 
numerically  more  efficient  than  the  conventional  methods 
that  work  directly  in  the  space  domain.  This  is  due  prima¬ 
rily  to  the  fact  that  the  process  of  Fourier  transformation 
of  the  coupled  integral  equations  in  the  space  domain  yields 
a  pair  of  algebraic  equations  in  the  spectral  domain  that 
are  relatively  easier  to  handle.  Another  important  advantage 
is  that  the  Green's  function  takes  a  much  simpler  form  in 
the  transform  domain,  as  compared  to  the  space  domain  where 
no  convenient  form  of  the  Green’s  function  is  known  to 
exist . 

B.  PURPOSE 

In  solving  the  open  boundary  slotline  problem,  only 
exponential  functions  arise  and  there  are  no  numerical  prob¬ 
lems  during  computation.  For  a  closed  boundary  structure, 
however,  hyperbolic  functions  are  required.  Knorr  and  Shayda 
[Ref.  3]  found  that  overflow  and  underflow  problems  resulted 
during  numerical  computations.  To  eliminate  these  problems, 
extensive  algebraic  manipulation  of  the  spectral  domain 
equations  is  required. 

This  thesis  presents  the  method  to  calculate  the  charac¬ 
teristic  impedance  of  fin- line  without  overflow  and  under¬ 
flow  in  equations.  Therefore,  this  thesis  is  a  direct 
extention  of  Knorr 's  work.  This  work  is  accomplished  by  lots 
of  algebraic  manipulations.  In  particular,  since  hyperbolic 
sine  and  cosine  functions  cause  overflow  and  underflow 
errors,  the  equations  developed  by  Knorr  need  to  be  put  in  a 
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=  o  (_o/j/y) 

£d  =  (PAXV^ 

,  /II 2JT  (jfc )  4>h  eUem 

.  (iiH)TT  (-§)  o44 

cv,d)x=  cv^pf  =  gujN-  (aftXWf-O  C%f 

(^p)i^(<4.o)1+  (*»f  DTAf-^J  CDA)a. 


C.  CHARACTERISTIC  IMPEDANCE 


The  definition  of  the  characteristic  impedance  for  an 
ideal  TEM  transmission  line  is  uniquely  given  by  static 
quantities.  Since  the  fin- line  supports  a  hybrid  mode,  no 
unique  definition  of  the  characteristic  impedance  can  be 
found . 

A  useful  definition,  however,  is 

-2-fav^  (eqn  2.70) 

where  V0  is  the  slot  voltage  defined  as 


2» 


( eqn  2.71) 
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where 

|  P  |  —  4n  (hi  -  On&n. 

=j  C^udxt  ~  • 


We  define 


a«=  (>,»+**(*«£*)] 

L  V/KHiX  /W34/J 

=  dn  cosk  £)i_P)  (eqn  2.62) 

1  '  V  *142  *(S^J  /*,!  J 

-  j  dii  oLp)  sroik  (?j9  (eqn  2.63) 


/«*43 


V/"U*  *3  4- 

=  ^Ui  cosU  ( &_P) 


/"U4 


)3 


(eqn  2.64) 


<U=  + 

V*46  ^53  ' 


~t  Max 


(/*g/*sq  _  W7  \  \ 

\/M6i  *?8  'J 


=  J  4ai  CViP)  50*1  h  QLd3 


( eqn  2.65) 


Also  we  define 
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(eqn  2.54) 


v  /»»+ 


v  -  ,*(53  <-  • 


Using  equations  (2.48)  and  (2.49) 


(eqn  2.55) 


(eqn  2.56) 


(eqn  2.57) 


Substituting  in  equations  (2.51)  and  (2.53),  we  obtain 


pc  _  7**9 

“  v*U4/*Ij4  ^46 '  U*U&*s*  / 


(eqn  2.58) 


Substituting  equations  (2.48),  (2.49),  (2.51),  and  (2.53)  in 
equations  (2.46)  and  (2.47)  we  obtain 


(eqn  2.60) 


(eqn  2.61) 


22 


The  dispersion  problem  is  now  reduced  to  the  form 


(eqn  2.45) 


where 


o^n  »d/jX 


A  numerical  search  for  the  value  of  which  satisfies  equa¬ 
tion  (2.45)  yields  the  propagation  constant  for  the  dominant 
fin- line  mode. 

From  the  equation  (2.39) 


ftflu  +  /Wu-  4-  =  o 

(eqn  2.46) 

A»tufle  +  m a3Ce  r 

/^x5  AK  4  Bh  +  M*\  CK  =  0 

(eqn  2.47) 

*nfi  fte  =  cf  S 

(eqn  2.48) 

m^i\  Pf  4-  =  if'Ejc 

(eqn  2.49) 

^33  Cc  4*  —  O 

(eqn  2.50) 

/m43  c€  -f  /m44  tf*  -+■  /nwu  -  O 

(eqn  2.51) 

A*»5t|Ch  4-  A*58  -  0 

(eqn  2.52) 

/v^^B€4-  /W*4-De  -yMtfch  -v  4»U8Dh  =  o 

(eqn  2.53) 

From  equations  (2.50)  and  (2.52) 
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Figure  2.1  Assumed  Electric  Field  Component  in  Slot 
in  x-direction  versus  x  for  Fin-line. 
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The  matrices  [He]  and  (HJJ  differ  in  only  the  last  two  rows. 
Each  is  a  square  8x8  matrix.  Using  equations  (2.39)  and 
(2.40),  we  may  write 


(eqn  2.41) 


From  equation  (2.41),  using  the  four  elements  in  the  lower 
right  hand  corner  of  the  matrix  MJ  ME  ,  we  obtain 

'  Exfc*)  3x6*0 

Ez(fa)  Tz&a)  .  (eqn  2.42) 


where  the  2  x  2  matrix  153  contains  the  Fourier  transforms 
of  the  components  of  the  dyadic  Green's  function  for  this 
structure . 

A  solution  to  equation  (2.42)  is  obtained  using  the 
Method  of  Moments  [Ref.  7].  For  this  problem,  we  have 
chosen  to  approximate  the  field  between  the  fins  as  shown  in 
Figure  2.1: 


e*00  = 


1*1  £x 


else-iOhc-^C. 


(eqn  2.43) 


e*(0  =  o 


(eqn  2.44) 
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eight  unknown  ;oef f icients  A  through  D  are  not  indepen¬ 
dent,  but  car  be  related  to  each  other  through  the  conti¬ 
nuity  conditions  of  the  field  components  at  the  interfaces 
between  the  three  spatial  regions. 

If  we  denote  the  Fourier  transforms  of  x-  and  z-directed 
current  density  and  electric  field  component  by 


=  F  { e*tt>) 

(eqn  2.35) 

EiU)  =  F  {  eaO0} 

(eqn  2.36) 

-  F  [  j„(*)} 

(eqn  2.37) 

(eqn  2.38) 

The  resulting  set  of  linear 

matrix  form  as  follows: 

equations  may  be  written 

(eqn  2.39) 


(eqn  2.40) 
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§*(*.,*)  -  B «U.Jyr4.Vii  +  C*CA)ajti  Vi1} 
ff  G/«,J)=  D‘&(«)5r*k  &Clu.+3) 

1^,3)  =  phwc^u  X  (oh, -3) 

$ z  {j^lP  —  Xjj +  t^C<6i)coslt  )ij 

$j  (^,3)  =  r^c^cMh 

where 


(eqn  2.33) 

(eqn  2.34) 

The  following  observation  about  the  solutions  in  region 
2  must  be  emphasized.  Any  wave  on  this  inhomogeneous  wave¬ 
guide  structure  will  partly  travel  through  air  and  partly 
through  the  dielectric  slab.  It  is  important  to  observe  at 
this  point  that  may  be  less  than  zero  in  any  of  the 

three  regions  of  the  structure  under  certain  conditions. 
When  o6n  =  0  and  k*  approaches  k0  (where  k«j  is  the  wave 
number  for  free  space),  ^  is  less  than  k^  and  so  Vf”  <  0 . 
Under  this  condition  the  hyperbolic  functions  in  all  three 
regions  are  replaced  by  trigonometric  functions.  If  k0  <  $  < 
then  vr  and  Vt  are  greater  than  zero  and  ^  <  0  for  some 
values  of  and  the  trigonometric  functions  replace  the 

hyperbolic  functions  in  the  spatial  region  2  only.  This 
suggests  that  the  nature  of  the  field  is  dependent  upon  the 
values  of  the  transform  variable,  .  For  the  conditions 
when  <  0,  yjj  replaces  such  that  () The 


(eqn  2.28) 
(eqn  2.29) 
(eqn  2.30) 
(eqn  2.31) 
(eqn  2.32) 


At  x  =  lb/2  : 

The  final  boundary  conditions  occur  at  x  =  b/2  where  the 

tangential  components  must  be  zero  in  all  regions. 

Ee;(lfc)M,*)=0  (eqn  2.22) 

=0  (eqn  2.23) 

B.  SPECTRAL  DOMAIN  APPROACH  TO  DISPERSION  CHARACTERISE 

The  scalar  potential  functions  can  be  transformed  into 
the  spectral  domain  via  Fourier  transform  defined  as 


(eqn  2.24) 

The  scalar  potential  functions  satisfy  Helmholz's  equations 
in  the  three  spatial  regions,  thus 

Vt  fa  t  ki  “  °  (eqn  2.25) 


where  denotes  the  two  dimensional  Laplacian  operator  in 
the  transverse  (x,y)  direction.  The  Helmholtz  equation  2.25 
are  transformed  into 


(eqn  2.26) 


where  =  aim  -  K*;,  =  »^n  +  p  -  K  .  Above  equation  has  solu¬ 
tions  after  applying  boundary  conditions  at  y  =  D  +  h|  ,  and 

y  =  -h^. 


J*  ^.,4)  =  (sW)  5T*U  V,  ) 


( eqn  2.27) 
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Similarly,  tangential  magnetic  fields  must  be  discontinuous 
by  corresponding  surface  current  densities. 


_  ( 1*00  eTi  l*Ui 


|Xl<x 


(eqn  2.14) 


1  o  |xK  X 


(eqn  2.15) 


AT  y  =  0  : 

The  tangential  field  components  at  the  interface  between 
region  2  and  3  must  also  be  continuous. 


Ew(o(,o,e)=  EasC*,®.*) 

C*>  o,  z)  =  H»  (*>  °> 

H -  HX3 


(eqn  2.16) 


(eqn  2.17) 


(eqn  2.18) 


(eqn  2.19) 


At  y  =  -t\2  ■ 

Once  again  at  the  shield  wall  in  region  3,  the  tangential 
field  components  must  be  zero. 


E?3C/)^  “}«.,£■)  —  O 


(eqn  2.20) 


£*3  COC,  ==■  O 


(eqn  2.21) 
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HJ= Cr-^  -)^^)ere 


(eqn  2.7) 


where  propagation  in  the  z  direction  is  assumed.  We  will 
also  assume  here  that  £,  =  £3  =  €0  and  £*  =  £®£r  . 

At  y  =  h|  +  D  : 

Applying  boundary  conditions  at  the  walls  in  region  1, 
tangential  field  components  must  be  zero. 


E»0X,h,+D,O=° 
E*|  (K,Ki  +  P;«)  =  o 


(eqn  2.8) 


(eqn  2.9) 


At  y  =  D  : 

At  the  interface  between  region  1  and  region  2,  tangential 
field  components  must  be  continuous. 


£ 21  (*  ,p  ,  2)  -  (*  ,  2) 

Exi  (1C,  P,  2)  =■  Exi  CK,  Pi  2) 


(eqn  2.10) 


(eqn  2.11) 


Also  the  electric  fields  at  y  =  D  will  exist  only  in  the 
slot . 


1*1 2  X 


/•  U  X 

=  j  . 

<  e,(oer 


(eqn  2.12) 
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o 

1*1 


(eqn  2.13) 
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II.  THEORETICAL  AT' * LYSIS  OF  FIN  LINE 


A.  FIELD  AND  BOUNDARY  CONDITIONS 


The  fin- line  supports  a  hybrid  field.  The  axial  compo- 
nants  of  TM  and  TE  modes  are  then 


E»  =  K?<t>e(*’P  er‘ 

Hi  =  kci<tT‘W.«erI 


( eqn  2.1) 
(eqn  2.2) 


where  the  scalar  potential  functions  ^  satisfy  the 
Helmholz  equation,  and  we  assume  lossless  propagation  so 
that  p  =  ±jp  . 

Further 


(eqn  2.3) 


with  kj,  -  ,  a.=  |  ,  Zt  3,  defining  spatial  region  of 
f inline  as  stated  in  Figure  1.1. 

Through  Maxwell's  curl  equations  the  transverse  field 
components  are  then  determined  by  these  axial  components  and 
can  be  given  as 


&=(r  'lor  efz 

Er(r*f +>>**£)  e>* 
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form  where  only  the  hyperbolic  tangent  function  appears.  The 
characteristic  impedance  is  computed  after  the  spectral 
domain  technique  to  find  the  dispersion  characteristic. 

Following  the  theoretical  analysis,  an  explanation  of 
the  computer  program  used  in  determination  of  characteristic 
impedance  is  presented.  Numerical  results  are  then  compared 
with  known  data  for  ridged  waveguide,  slab  loaded  waveguide, 
slot  line  and  fin-line  [Ref.  3].  This  comparison  estab¬ 
lishes  the  accuracy  of  the  numerical  results. 


w  - *r  ’  V"’7>T.:  VT 


ex(x)  is  arbitrarily  selected  as  1/W  so  that  W*ex(x)  =  1. 
PAV^  is  given  by 


EX h*  ■  Si  A(\ 


S  ■ 


(eqn  2.72) 


Pasval's  theorem  is  applied  to  eq(2.88)  to  obtain 

4^. 


This  expression  must  be  evaluated  in  each  of  the  three 
regions  of  the  fin- line  shown  in  Figure  1.1  Therefore  the 
power  flow  may  be  expressed  by 


i  Jr.^CP'  +  P* +  PO 


(eqn  2.73) 


Since  P*<J  can  be  determined  after  finding  the  value  of 
^/A  •  Equations  including  only  hyperbolic  tangent  functions 
and  the  slot  voltage  V0  have  been  developed  the  character¬ 
istic  impedance.  The  lengthy  algebraic  manipulations  are 
shown  in  appendix  B. 


I 
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III.  COMPUTER  PROGRAMMING 


A.  NUMERICAL  ANALYSIS 

The  computation  of  characteristic  impedance  is  based 
upon  the  solution  to  the  dispersion  chracteristic  problem 
for  the  fin- line  under  consideration.  In  other  words,  the 
wave  propagation  constant  or  wavelength  ratio  X  must  be 
known  before  any  other  investigations  can  be  started  since 
only  in  this  case  are  the  scalar  potential  functions  in  the 
transform  domain  known.  The  computer  program  for  the  wavel¬ 
ength  ratio  X l\  is  already  developed  by  Prof.  Knorr.  The 
next  task  is  the  preparation  of  the  appropriate  equations 
for  the  numerical  evaluation  of  the  time  average  power  flow. 
For  ease  in  numerical  calculations  and  for  programming 
purposes  all  geometric  parameters  are  normalized  as  follows; 
hj/D  ;  fin  location  relative  to  the  positive  "y"  side 
wall  normalized  with  respect  to  D 
hz/D  ;  fin  location  relative  to  the  negative  "y"  side 
wall  normalized  with  respect  to  D 
b/D  ;  waveguide  height  normalized  with  respect  to  D. 

It  is  observed  that  for  the  power  flow  computations  the 
infinite  integration  in  equation  (2.88)  is  replaced  by  an 
infinite  series  due  to  the  discrete  nature  of  the  transform 
variable  a (/*  .  Preliminary  numerical  investigations  of  the 

coefficients  of  these  series  indicated  an  even  distribution 
with  respect  to  the  variable  so  that  computation  time  can 
be  saved  by  summing  over  a  half  interval  only. 

The  time  average  power  flow  equations  are  prepared  for 
six  different  cases  as  follows; 

case  1  ;  (V,Df  <  0  in  region  1 

case  2  ;  (ViD)2”  >  0  in  region  1 

case  3  ;  (Va.  T>t  <  0  in  region  2  and  (V,D  f  <  0 
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case 


in  region 

3 

4 

;  (VxDf 

< 

0 

in  region  2 

and 

C&D  f  > 

0 

in  region 

3 

5 

;  ()Uf 

> 

0 

in  region  2 

and 

(V»D  f  < 

0 

in  region 

3 

6 

;  (^Df* 

> 

0 

in  region  2 

and 

(Ad  f  > 

0 

in  region  3 

After  this  preparation  of  the  equations  for  the  computer 
programming,  the  only  remaining  task,  is  to  investigate  the 
numerical  .  integration  with  regard  to  its  limits  and  the 
procedure  to  be  used.  One  of  the  limits  of  the  summation 
discussed  previously  extends  to  infinity  so  it  is  necessary 
to  determine  an  appropriate  at  which  to  truncate  the 
computations . 

From  the  representative  examples  shown  in  Figure  3.1,  it 
is  seen  that  the  coefficients  in  this  infinite  series  are 
found  to  decay  rapidly  so  that  a  finite  approximation  yields 
good  results.  Figure  3.1  shows  the  characteristic  impedance 
as  a  function  of  n,  the  number  of  terms  in  the  truncated 
series,  for  various  slot  widths.  The  value  n  =  50  is 

sufficiently  large  to  be  considered  infinite  for  all  prac¬ 
tical  purposes.  It  ’  *  noted  that  when  w/b  =  1,  Ex(«6n)  =  0 

for  n  >  0.  Therefore  whenever  w/b  =  1  only  the  n=  0  term  is 
computed  in  impedance  calculations. 

B.  COMPUTER  PROGRAMMING  AND  LIMITATION 

The  first  part  of  the  FINIMP  computer  program  finds  the 

/ 

\/\  which  makes  equation  (2.45)  equal  to  zero.  The  assump¬ 
tion  is  that  A/ A  is  between  0.1  and  3  in  the  giga  hertz 
range  and  continuously  decreases  as  the  frequecy  is 
increased.  Since  this  program  is  very  sensitive  there  is  a 
possibility  that  the  wrong  value  is  searched.  In  that  case, 
we  can  easily  notice  because  there  is  abruptly  jumped  value. 


At  that  time,  if  the  assumed  range  of  X/l\  is  resticted  to  a 
smaller  interval  than  before  the  correct  value  can  be 
obtained. 

The  rest  of  the  FINIMP  computer  program  is  the  computa¬ 
tion  of  characteristic  impedance.  Since  all  equations  are 
composed  of  hyperbolic  tangent  functions,  there  is  no  over¬ 
flow  in  the  computer  program.  But  if  (Va  Df  <  0  tangent 
functions  replace  the  hyperbolic  tangent  functions.  At  that 
time,  there  is  a  possibity  of  overflow  because  tan  nir/2  is 

infinite  for  n  =  1,  3,  5,  .  In  that  case  the  tangent 

function  should  be  set  to  some  value  which  is  the  maximum 
value  of  computer  ability  (46**  (1-16”®)  >  |tan(x)|). 


IV.  NUMERICAL  RESULTS  AND  COMPARISONS 


To  check  the  accuracy  of  the  numerical  results  gererated 
by  the  computer  program,  comparisons  are  made  with  data 
available  in  the  literature  for  ridged  waveguide,  slab 
loaded  waveguide,  slotline,  empty  waveguide  and  fin- line  as 
outlined  in  [Ref.  3]. 

A.  RIDGED  WAVEGUIDE 

When  =  1  and  w/D  <  1  or  when  the  dielectric  substrate 
thickness  D  is  reduced  to  zero,  the  fin- line  structure 
becomes  ridged  waveguide  with  zero  thickness  ridges. 

The  impedance  of  the  ridged  waveguide  has  the  following 
relations  [Ref.  8], 


2?o6o 

0  -  cfafj*- 


(eqn  4.1) 


When  the  guide  width  and  the  slot  width  are  equal,  the 
ridged  waveguide  becomes  an  ordinary  rectangular  waveguide 
for  which 


AXc 


\l 


2JX 

b 


(eqn  4.2) 


When  the  guide  width  is  increased  at  fixed  slot  width,  the 
characteristic  impedance  at  infinite  frequency  and  the  free 
space  cutoff  wavelength  go  asymptotically  to 
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(eqn  4.3) 


The  cutoff  wavelength  ^  was  determined  from  [Ref.  11]  for 
several  values  of  W/b. 

Equation  (4.1)  was  then  used  to  calculate  the  character¬ 
istic  impedances  and  these  values  were  compared  with  the 
results  using  spectral  domain  method  of  the  computer  program 
discussed  in  chapter  JIL  -  Two  results  are  very  close 
together.  Figure  4.1  and  4.2  shows  the  computer  results. 
These  figures  agree  with  the  previous  Knorr  and  Shayda's 
results . 

B.  DIELECTRIC  SLAB  LOADED  WAVEGUIDE 

A  variation  of  the  fin- line  structure  of  Figure  1  where 
w/b  =  1  and  £r  >  I  results  in  a  dielectric  slab  loaded 
rectangular  waveguide.  The  slab  loaded  waveguide  has  also 
been  studied  by  Vartanian  et  al  [Ref.  9] .  They  consider  a 
guide  with  the  slab  centered  in  the  waveguide  and  they 
present  an  analytical  expression  for  the  voltage  impedance 
Zpv  at  the  center  of  th -  slab.  The  impedance  computed  in 
this  analysis  was  spe<  iied  at  the  edge  of  the  dielectric 
slab.  Therefore  a  relationship  between  the  two  impedances 
must  be  defined  before  any  comparison  can  be  made.  The  impe¬ 
dance  at  the  edge  of  the  slab  can  be  easily  defined  by  the 
expression 


(eqn  4.4) 
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SLAB  LOADED  WAVEGUIDE 


© 


Figure  4.3  Characteristic  Impedance  Z  vs. 
Frequency  for  a  Slab  Loaded  Waveguide. 


35 


where  E®*JC  is  the  field  at  the  edge  of  the  slab  and 
is  the  field  at  the  center  of  slab,  thus 


(eqn  4.5) 


where  the  various  quantities  are  defined  in  [Ref.  9]  as 


5  = 


(eqn  4.6) 


(s')  -  £r(1v>)  "(lij)  (eqn  4.7) 

<f-  =  (aT)i('*)  [^V  ~  (  «T}  (eqn  4.8) 

and  c  =  D  is  the  dielectric  slab  thickness. 

The  characteristic  impedance  of  a  slab  loaded  guide  is 
computed  using  the  spectral  domain  method.  The  results  were 
compared  with  equation  (4.5)  results  and  [Ref.  3].  Good 
agreements  are  obtained.  Figure  4.3  shows  the  slab  loaded 
wave  guide  impedance  . 

C.  SLOTLINE 

If  w/D  <  2  and  is  sufficiently  high  for  the  fin-line 

structure  of  Figure  1.1,  the  field  is  tightly  bound  to  the 
slot.  For  this  condition  the  presence  of  the  shield  will 
have  little  effect  if  the  walls  are  sufficiently  far  removed 
from  the  slot.  In  this  case  the  fin-line  will  behave  like  a 
slotline.  This  behavior  is  illustrated  in  Figure  4.4  where 


37 


the  characteristic  impedance  of  a  fin- line  with  w/D  =  1,  £r 
=  20  have  been  plotted.  Also  these  results  are  in  good 
agreements  with  [Ref.  3]. 

D.  SEVERAL  FIN  LINE  IMPEDANCE  CURVES 

Fin-lines  are  generally  enclosed  with  a  shield  that  is 
compatible  with  the  dimensions  of  the  standard  rectangular 
waveguides  for  the  millimeter  wavebands.  Above  22GHZ,  all 
these  guides  have  aspect  ratios  b/a  =  0.5.  Further,  the 
fins  are  most  often  centered  in  the  guide  and  printed  using 
D  =  0.005  inch  substrate  with  =  2.2. 

Figures  4.5  -  4.12  show  the  impedance  computed  for  fin- 
line  v ith  WR(28),  WR(19)  and  WR(12)  rectangular  waveguide 
shields.  These  figures  may  be  compared  with  data  for  the 
same  structures  as  presented  in  [Ref.  3].  Such  a  comparison 
shows  excellent  agreement  for  all  but  the  small  values  of 
W/b.  For  small  values  of  W/b  (W/b  =  0.1,  W/b  =  0.05,  W/b  = 
0.02),  the  results  of  the  FINIMP  program  show  lower  impe¬ 
dances  than  the  results  presented  in  [Ref.  3].  As  the  W/b 
is  smaller,  the  difference  is  larger;  for  W/b  =0.1,  2-5 
ohm  is  lower,  for  W/b  =  0.05,  9  -  15  ohm  is  lower  and  for 
W/b  =  0.02,  19  -  25  ohm  is  lower. 
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Fieure  4.5  Charact e r 1 s t 1 c  Impedance  Z  vs.  Frequency 
for  a  Fin-line  With  b  D-28.0  h|, D-28.0  h,£/D=2/.0 

D-  .  005 "  y-  =2.2. 
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The  matrix  elements  of  ("Mg)  and  fMj}  are  normalized  at 
this  point  with  respect  to  D,  the  dielectric  substrate 
thickness.  The  normalized  matrix  elements  are  now  presented 
in  two  forms.  The  first  of  the  element  equations  is  for 
(Y;  Df-  >  0  and  the  second  is  for  (ViD)2"  <  0.  For  the 
matrices  fMg]  and  C^j3  »  the  elements  m||  through  m^8  are  the 
same . 
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APPENDIX  A 

SPECTRAL  DOMAIN  MATRICES 


The  continuity  conditions  are  transformed  via  equation 
(2.24)  into  the  two  dimesional  Fourier  domain.  The  solutions 
to  the  two  Helmholtz  equations  given  by  equations 
(2 . 27 ) - (2 . 32 )  are  substituted.  Finally  a  matrix  form  of 
linear  equations  is  derived  as  follows; 
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B.  RECOMMENDATIONS 


Coupled  fin- line  will  find  future  use  in  building  direc¬ 
tional  couplers  and  filters.  For  this  purpose,  the  normal 
mode  wavelengths  and  impedances  are  required.  A  special  but 
important  case  is  that  of  symmetrical  lines.  Coupled  fin- 
lines  for  which  the  normal  modes  are  odd  and  even.  The 
program  described  here  may  be  extended  to  cover  this  case  by 
following  the  procedures  outlined  by  Knorr  and  Kuchler  for 
coupled  slotlines  [Ref.  5].  This  should  be  accomplished  to 
improve  the  utility  of  the  program  described  in  this  thesis. 


V.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  spectral  domain  technique  used  in  conjunction  with 
Galerkin's  method  has  been  presented  to  calculate  the  char¬ 
acteristic  impedance  for  the  dominant  mode  of  the  fin- line. 
It  has  been  shown  that  a  matrix  formulation  of  the  problem 
permits  the  elements  of  the  dyadic  Green's  function  to  be 
calculated.  Solving  the  matrix  equations  leads  to 
containing  only  hyperbolic  tangent  function.  These  equa¬ 
tions  circumvent  the  overflow  and  underflow  problems  which 
occurred  using  the  previous  formulation  presented  by  Knorr 
and  Shayda. 

Numerical  results  obtained  using  this  method  have  been 
presented  and  compared  to  other  existing  data.  Good  agree¬ 
ment  has  been  obtained  in  all  cases  thus  establishing  the 
accuracy  and  applicability  of  the  method  for  the  full  range 
of  structure  parameters. 

There  is  a  possibility  that  tangent  functions  may  cause 
an  overflow  problem  if  (V;  0  'f'  is  less  than  zero.  In  this 
case,  the  value  of  tangent  function  can  be  obtained  within 
the  capability  of  the  available  IBM  3033  used. 

In  this  thesis,  paticular  interest  is  devoted  to  the 
computation  of  fin- line  impedance.  Fin- line  may  exhibit  the 
characteristics  of  ridged  waveguide,  dielectric  slab  loaded 
waveguide,  slot  lines,  and  conventional  rectangular  wave¬ 
guide.  All  of  these  structures  are  fin-line  substructures. 
So,  the  computation  of  fin- line  impedance  permits  all  of 
these  structures  to  be  analyzed. 
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Figure  4.12  Characteristic  Impedance  Z 
for  a  Fin-line  With  b/p=18.8  h|/D=28, 

D=  .005"  =2.2. 


vs .  Frequency 
2  h*/D=8?4 
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FIN-LINE  WITH  WR(12)  SHIELD 


IMPEDANCE(OHM) 

JOO  120  140  160  180  200  220  240  260  280  300 


FIN-LINE  WITH  WR(19)  SHIELD 


Figure  4.8  Characteristic  Impedance  Z  vs.  Frequency 
for  a  Fin-line  With  b/D=18.8  h|/D=18.8  h2/D=17.8 

D= .005*  tr  =2.2. 


FIN-LINE  WITH  WR(28)  SHIELD 


FREQUENCY  (GHZ) 


Figure  4.6  Characteristic  Impedance  Z  vs.  Frequency 
for  a  Fin-line  With  b/p=28.u  h|/D-28.0  hz/D=27.0 

D=  .005  £r  =2.2. 


w£  =  f  ««PJ>CfP)5r<*ih  rev,  DXhX)i! 
i  j  C^0)(pDV  srmCcVeDX^J 


■W/j2  —  —  Q 


f  r  j(^yWP)CV,p)5whCCViP)C%)3 

/Wng  =  < 

t  -j  ((4«DV  cVfD)  srM  CcV,"PXh/o)J 


^74  —  /W17  y  =  =  O 


£  _  f  Oi.Py's(>kCOlPX,’^3 
'Wil  I  rctf'px,%)j 


€_  £  ,£  £.  f  C  e 

/W&z  =  /W$3  -  /t>1f4  =  /W$£  =  M&b  =/vA&%  -  o 


T  .  y  T  T 

—  /Wxjz.  =  /»^y3  =  AtArfa  —  O 


x  f  C(4,pXc«hCcV,p)Cl’/pV] 


J 

AA\yS  - 

V 


Ci4ip)ic»src^pve|'^>)] 


54 


OiD) 

■  -jCfaof  srm  (Y/tO 

;  -  (&j.of-c°sU  cLs>) 


-j  (u)£,  p)(y,D)  cosh  CcV,DXh/D^3 

■  (to£,D)cf(,p)  c.cs  ccVi*»>x|yp)3 
-j  (iO£iP)(ViP)  c»shC&P) 

o o£1i>)alrp}e°«aL'b) 

-j^fiPXtosr^UC^P) 

j  (iO^djcV/d)  sr/*»08ft0 

5 

Ccy, dXVp)3 
.  (o&Dy^D)  cos  ccVfDxypX) 


APPENDIX  B 


& 


K  ■ 


( 


TIME  AVERAGE  POWER  FLOW 


The  following  expression  for  the  coefficients  A^  through 
are  derived  from  equations  (2.54)  -  (2.61). 

If  (Y,  D)1  <  0 


"J  ~ 


c  (V,#dXk/d)J 


] 


,h  =  rr  CfriPf P*~Ex  -6kPXftD)p%  1 

J  L(w/ap)ckci  d^d)  s  J 


If  (  Y,  Df-  >0 

aC_  _ 

Ah=  -  i  f  -  G^PXfoOCfeg  1 

^  L  (^PXV.oXfoeX^h  CcY,PX%)J J 


If  OLDf  <  0 


r®’—  ~  dzxD  £& 

dei  cosCVaP) 

A\  >  f  oLi  -  d»  t^Ex  1 
j  L  det  CY/p)  Spn  J 
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dz\  d\Z 


where  det  =  d,,  d^^ 
If  (V^D)2-  >  0 


(f  =■ 


*•  ■}  G 


du-D  Ex  “tUi-D  Ez 
c$€t  cosh  CVaD) 

%  r  d«  "cU  *| 

J  L  cjct  &(Os6nk(to)j 

<  0  and  (VjDf*<  0 

.  f  0>£*D>CV/d 

1  rdu 

loa^DXV'/p^oy^ ro^b)C*%)J J L 

.  r&MDYpvCCKc^-terffl  1 

kde*OW))sc*»i)ifb)^ 

GUDXpp)c<:&*D^M>i]i  r 

<JuI?£c-cLl&'| 

.  Ck)/*D)CY/P)(Mt>  f-  Jl 

jdcccsCKi b)  J 

***  ro&x %)3i rd,,  d&  -  ifg  i 

CWdcjPf  JL 

dctOL^rm^J 

>  r  ^  -  daltfSi')  "I 

ckt  CCS  OSTp)C MO**1* CCVrt>X%>J J 


J 


d<£  C  HMsr^oLrPX^0f-cos^  £0>*t 


If  (YxDf  >  0  and  (Y3Df  <  0 


f  Ole&KK pK£J>?  1 

“<JU[)£a  1 

m>»ob)<&>otjb«\  [cYj^x^yiJ 

1  dttcashQLD)  J 

.  f  ]  r  JnC&  -<Lds k  i 

Tl  (^vOLPXfaPf  JL^tO^>)spvih(toJ 


„h  [GmV)(^  fO^#-  0<sDfJ  -i  pejp.cfS-cU^I 

p  ctJ/«3){%.p)^PY  JlJ*  C£?3  KOip)  J 


r  (Ys  p)(KcxD^tv»\  rcy/px^)3 1 

L  ciJooteDf  ) 

nf  =  -  f  G^DfCduCfg^-jal 

^  Lde> c/Jo Q^ofsxl*ir(yfp{(jy£)')} 


pH  _  '  r _ c^zp^cincfgK  rf£j') _ 

j  L^eta^srmH  aj>x^DY"cDs[;(y^p)c»%)3 


If  (  YaDf  >  0  and  ( Y3Df>  0 
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_  r _ QJt&'xSiOXfavf' _ *ir  ducfex -<k±&£z‘ 

D  lC«ed>;Oj))(XcsD/*^h C&WCWjJl  detcosKOU))  _ 


+ 


7^dX(sp) 1  f  ^ntf^-^pS  1 
.  (U6*D)&DXfoD/~  JLaetOlP)Sw*»hC^DXI 


J 


"T&DVgp)  fflbfif-fepfi  1  f  dbjf£x-tiaCf&  1 

.  J*'  de*e*sKO^»  J 


c^OLp)s.>K(y^ 


(^f(duQx&-(Ui?&') 

detc»sUQ/zVXp;*0fSV\\*  CViDX^J 


■NT 


CKo-Pf'CAt  C^Ex  ~  Am  D^Eg) 


de*  Sf/nK  0,xP)(X5D)xccsh  CO^X^J 


a 


From  the  above  equations 

if  (V,  vf  <  o 

gift  _  f;T 1 1 W[CK,pfC%^1 

D4  (fei  q)4  favfl1[cy,-,D)C%)3 


61 


i #  (fit.  iteiPfa  -a,PXW&f[it**nxCtfc>xl%>)l 
I +  ■  (uWof  Ol'pfOi.DJ4  ^  CoV'dxVd)] 


/r-rtf _  r,  ctfapfe  -6hcocgio&Cn-*>n  Cb|.'i>xl^a3 

D4  (iJ/lP)(V'D)(^aD)rt 


(iff  Ah  _  gjjjf 
D4  C4 


If  ( V,  D J2-  >  0 


Ae-ifff .  f?  ci -jba^rtpx^ai 

D4  05uD)4*wiW1tcy1PX,^fe)3 


_  feuPfS  -  ^p)gp)&ft-Awi'Jrol,pxlix>J) 

pa  C‘JA«pf'apfoiio)4*»«hicci(.px|'>6a 

^  tff  .  fa  p)^  -tt.p)(8p)MP -wrtrtpxtea] 
p+  J  L  Cfe.p;4(>*'/«p)CV,p;*«ii,ifcv,p)c%ii 


art  p£  _  j£j£f 
D4  p4 
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_  rl 
D4 


=  Hk*  J  [i+W-oSfpjJ 


cj^f_  rJ.S-cLS 
D*  l  deiOt’c)  J  L  WO/P)  J 


Ccl‘)^_  .  f'CdaS-Ja&')^.r«-iig)~irH-W(Vxfe)  1 
D4  J  L  dei1  Jloi$>JMGtf» 


cK  _  _  c«  (c^ 
D4  P4 


(ViDf  >  0 


Ch  (cfcf  _  rdn  £  -  da  fe  T*T  I- jbuiK^CVj))  ~i 

D4  L  drtclW  J  l  jbwl^cVj))  J 


tlfch  _  c Hit 


If  (Va-Df-  <  0  and  (V3Df*<  0 


r  rtJ£,pyy/p~Kfcj>f _ frc^i 

p*  U‘J^)0'PXfcoWoj*M<[irv3'DX^]J  l  d+  J 


Cue»ox:Vj,py(fapf 
ciJ^^oa'px^iP/cvA)  Jwi  Ccv/ p>?% j 


r ouwwor rfcof-  i  r  (*  (<*fi 

l  CMit>XW(faof-  Jl  p*  J 


,  r  ffop) -(fe^pfj  fp*  fc^f i 

L  ((o^^cfepy  Jl  p+  J 


&_(££  _  r  (4p)fftp)r(t^pf-cfe.pf]-i2-[-  cg  (yy~i 

P+  l  C-WOiftOlJoD)1  J  l  p+  J 


,  r  66,pygp)  Rfapf-  (fe»pf3  ~|  . 

4 1  c^Dxt'Px^pr  J 
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-ir(cty-.c>'  i 

L  ClSoXVaDf  JL  d4  J 


bc  (&?  _  .  r  ciXiP)(ypx:faD7>~ _ 

p4-  J  lClO£^(i'/C))()iiDX-^&V/D)Ch?fe')J 


~^pXfapfi»mrol>pXl%0 1  r  ch- (ty i 

o£px*»p/  JL  p+  J 
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(&f-  &  _  &-  (&? 

P*  "  {A 


gl  (£T  r  c^pYiMcKoPf  irce(ckf-i 

p+  '  uooilwawtfoo^Ewc^jL  jD+  J 


r^D)ftp)r(^f-Cfe>Dt3]r  d>  (c^fi 

L  3L  p4-  J 


(gy  t?  (cyf 
D+  P* 


ef-gy  _  .  r  rw£jP)Q,irp)ct^-Df'  ir ^  iff] 

p+  Jl^^xv/p)fepfi«»>cof/ox%)JJL  p*  J 


J  l  oi^ai.-'pxfooy  JL  J 


(fXc*  _  _  &-(c*f 
a*  pf 


&-cef  _  r  Gt»py»p)Cirfopf-  4apf3 1  ~i 

p+  L  cu/»p)(&'pxn»py  JL  p4-  J 
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,  r cVibyfapfjwDnbbx^i  i [■<*■■  o*f -1 

1  OSTcXISaP)1-  JL  jp4  J 


<$f£  = 

P4  P4 


_  -  r6t«PXfrPX0fap/-(faPfJ  'lf'd^g'1 

D4  J  L  (w/tPXY/DXfoP)1  Jl  j  D4  J 


, .  rcv^cfapf^  rAtoc%^i  rch  i 
JL  cVz’px.fc>py-  JL  d4  J 


if  (tf)*  r  (fapf  fr  t+ Wfdi.V%>]  1  r  tf  iff  1 

pa  Lcfop)*  J  L  WRy/pxVpy  J  L  d4  J 


1^£= [«§>’  ^Kfecyp>]j[-^j 


£jjt  r  ffapf  fr  i+ W-CcViPX^I  -I  r  ce  (c^i 

[)4  L(foP)1  J  L  iamCO^X^J  JL  JD*  J 
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For  the  case  (VjD)2"  <  0,  V,  is  imaginary  in  which  case 
(Vj^D)3"  =  -(V,  D)z.  The  power  flow,  P(^  ,  for  this  case  is 


CfD)(UAD) 


ceo)&j&  -  Qpj«.o)^gfie  -  (wfouv 
[*<*'«(■£) -sw»  [^Ctt»4)]] 


+ c^)m,o)d’o)  [LM.  t0c,Df- 

«.»)  (J^pr  ]  [  2  oToc1)*) + s™ 


REGION  2 

For  region  2,  the  power  flow  expressions  are  the  same  as 
for  unshielded  slotline  except  that  the  Fourier  integral  is 
replaced  by  a  summation  and  the  interval  2  is  replaced  by 
the  interval  b. 

For  the  case  (yz  D)1  >  0,  Vj.  is  real  and 
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D+  D+ 

REGION  I 

The  term  P|A  will  be  used  for  the  power  flow  in  region  1 
for  the  case  (Y.D)2-  >  0. 

pi  a  =  -  i(i)Pe  £"  \  UpKMWffif 

_  j  (pf^pxy.p)  +  apf«.woiD) 

Atf-Ahr  ,/h.vi  i  Om 

c*  J1  aoTC-wrrtox®  2(d)JHQ« 
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W f  J"  ■'T!' W TT 5-T  ^ T wv  ■’  w -;  ^  •—  TT V"  Cl °2T  H‘"  V' "Tr-T 


,  .  rctiPXteipfjbwt'  ra>p)(t%a  ir  cty  c1*  ~ 
rjl  CtPXI4»BjL  JL  j'0+  . 


<*  _  _  ^  ■  C<*t 

p+  p+ 


e>*  (ch)^  pa.oV8P)KK3Pf-  Cfe?pj3  ~i  nyfcK  i 

p+  l  cu/wo)OtD3045Dy-  JLjd+J 


,  rd>p)<,fej)f^wt' CcV^-fl 

l  OPXfoof  JL  P+  J 


(E,l'fe  ch  _  Q-rf 

D+  "  p+ 


pe  ■  (ft  _  roapff-f ]f~cg  fc*f 

P+  L  CKiPf -1  l  jfeM.hlCcY,»(l%)}JL  i*  . 


£i^f  _  -  rffaptf ri-^^cto66jirce  (cfcf i 

0+  HoWfJ  l  JtrvnU & v„o > J L  jp+  J 
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&■(<?■?  _  -  r^ov^nfeDf-oao^irce  (^i 

p+  J  l  c^px^xti^y-  JL  p+  J 
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P)  (Kgj/iniih  fc 

OfiP)Oi3.D>2- 


,  r^pX^Dftbv^k  rokvtH))]  i2 rch  (C.hfl 

L  J  l  o*  J 


B*-(BKf_  >  ro^pY&pfO^Pf _ ; 

D+  J  L  cneoPXYj>)(feDji:"3bM.h[(;V>DX,%)J 


i&PWpjgl^f-fopfl-i  r<*  fcef -[ 

(iJ^UOiClPKHiD)2-  Jl  p+  J 


.  •  r«.pf( ttjfTCfaD^OoDf Jl  r C€  <<*)*  -] 

J  L  c^iPX^wokDXo^DX  Jl  j  p+  J 


-  rO^.PYfeP)ttfaPf-Oi:?Pt3i 

Jl  (u^XlraCfePf  J  ' 
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\  I  *V  •  ^1 '  V  F  ■  V  ^  ^  ^  ^  v  t  ,  —  \  w  >i  *■  -'  w 


t*  _  D tjgf 

p+  p4 


If  (ViDf*  >  0  and  (%T>f>  0 


£■■(#?_  | 

r  ci^iPxyjWcuicf-  f 

fa-CO*! 

Dt 

UuXiPXyj’Xt4»D7*w'KCCW>X!%)]J 

l  D*  J 

+  2 


"  fujgOOCfepXkiPf- _ 

-6J£2p)cvj>)0i?Df*».K[]:y,Px,%ii 


.  rc^mWc^f-C^tfvr d*  (frjh 

T  L  (.UliPXtPXk#?-  JL  [4  J 

b1*  {££  _  ri^PXf»p)SiiPf-<fe»pti  V-fcj, c^fi 

t*  *-  owdkPXMj*-  J  I  p*-  J 


.  r^PXWDfcJf-O^Pt) 
L  (u)W(Vj>XK.jp/' 
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+  -  rcv^pxfeipfja^rcv^)^)]  lr^  6>y-i 

J  ^  CY1PXfoP>3'  -it  D+  J 


(&f  c£  _  -  6*  •  c*)* 

P+  p+ 


ft  (cKy  rr^px^pjfo^pf-  fepf]  irctftr  eft] 
p+  t  (u^Dxfroo^py1'  jp+  J 

_  r  rt/pXfci#*™  KVfox:1#)]  1  fch 
•-  OPXftbP)*-  1  p*  J 


gs 

P* 


~irce  (t^fi 


.p^epXcfapf-O^pthrc1'  (c^i 

c^xiwcKcjpy-  JL  p+  J 

&■(<*)*- 

D+ 


(iztioxiioy&ipf  ir<* 
Cu)g1P)0tPXti3l>Xjhw>Q^,PXhjt))jll  p+ 


u.p'ifeff)  rc^xt>f- 1  fW  cjM 

CiOSiWOLDXfe^P)2'  Jl  jp*  J 
b *J£t 

pt 


64M>XtoXfc>P?-  Jl  D*  J 
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roi»ckj>ffa«  r 1 

L  (iLoXMP)’-  A  p4  J 

.  rOl>,PXI<-xDf*wiCnf3<bvi%-)]-|J-r<->i  (t^fl 

*-  otwcfeD)3-  J  L  J 

ii (£?-  -r  -  (^vK'DY^Dt  -. 

p^-  J  ccy»,D)c*^i>2l  J 

Rfap/-<fc<p)*-1 1  rye 

L  <r^p)oiD)cfejPX  Jl  p4-  J 


+J'RS 


PXVa-P)  (KuO~)4 
DXV^Oi^D^ 


J  L  fU)£iPJ(i0/l»D)OllP7,(fc1D)4  J I  p4  J 

-.-■  r  C4.p)Cgp)Rfapf-<fac)^  -[ 

J I  ('OfxP)OiP)tfoP)2‘  J  ’ 

rcV/DX^iP  r  <* .  (c^i 

l  cit-xKijpy-  31  p*-  j 


m 

V  v 


& 


I 


(g*f  Dh  _  [^(C^f 

”  tn 


If  (Va-D)3-  >  0  and  (Y3Df"<  0 


&’■(&)*  _  r  6&>t>XV>Wfcipf _ frc'CcTi 

P+  L  o3&waEiXfe^^Rto?%aJ  L  c*  J 


.  ,  r  fdfopxy/px^f- 
T  4  Lc^PXlrxfoo^jUn  COStoOWJ 


rcumoir rfojf-ffopfJ  ir^y*  <*  i 
l  WkeXlPXfcjDJ3-  J*-  c*  J 


+  r  a.g)(e>p)Co^pf-(fe^f3  -i^rcV^fi 

TL  (i0ej>)(tp)OiiDf-  J  L  p+  J 

bVb|^_  r^pygp^CO^pf-Cfap/J  i^-rcf^i 

p+  "l  Cn)«Dxyi.p)((t5P)  J  L  p+  J 


, :  r(^DX&t>)r<3^pf-o^ipfJ  ] 

l  6JJUD)CkD)Qi}D')  J  ’ 
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.  •  . ’  - '  - 


(£ f  B*1  _ 
I>+ 


B*  (£t  _ 
D+ 


P+ 


iLSf - 

D+ 


(E*f-Ce  _ 
t* 


P+ 


D+ 


C^PXttPXfciDf  TT^fc^l 

oj^ppcVxDjcfoD^jhwUCajKX^)^  i  d+  J 


rr^DYgpi  r(b:>pf-(fcDt1  irch  fch7h 
L  JL  qA  J 


pe.  (<*>)*■ 

p+ 


-  p  ('i04p)(V^px:fe»pyL 

J  U^pxy^xm>7*>*1>  Cd.px'tf#  p4-  . 


-  r  a.pxp«  r (fctpf-fept]  ifty  -  ch  i 
j  l  (kXiDXTOCteP)1  JL  J  [A  J 


^  (c8)* 
P* 


r(4,pxgp>  rcfapf-  qtaPf J  i  r  cc  •  <yy~i 

L  Cd/PX&'PXM)*'  jL  c*  J 


70 


amniaflaan 


u)£3d)CV3DX^d 


oj^p)  cTOc^py-Jkwc^Ji  jp4- 


&D)tf«Ca^f-(MLJ 


O^DXWCfc*^ 


LimaPtMi  g^px^)3  i  r^  &T 


CV/pXfoDr 


- ,  '•  \h 


nl 


_ CulSaOXVjDV^/- 


Cw€^X)i1))Oi»Ki‘AwihCDiPX>7'p)J 


l! 


cMrtxifoxk&f 

&  (£t_  rfaDXto'Xttatf'-  (Mft 


^mcfftook^r 

zGUOXW)  CC  fopf-  C  tevf} 


G*MO)CY£pX  favr 


=  -  8  (t)  R<  1  -Ljpp  +  CfPX^P) 


+  (^)C*JAt>) 


+  QOCrtiDXl'W  ^0-  +QPXWAliP) 

<«»  c±f-  +  (f*w»  [*&-  + 

Q&tCJmti) f ^-h  - 


I'WXfepK./.P!1  c*-(f*r  &  t<yfm 

l  CYi(»  D4  t  C4>*<0  D4  + 


(?#«)(»)  +  (^D)(um2Dji(y/D)  4 

-  £f*J 

f*  CM»  +  yr<«  (DSfofl  +iC(pX^W 

lowttwpsti-  «££][, -«<»*,>]} 


REGION  3 

For  the  case  CfcDf  >  0, 


is  real  and  the  power  flow, 


(Y'd)  Jt$£  +  (f  »^w  + &&**>*£&) 


pb  (W 


jpr  +<swe»x®)&££-  -(fpfu.»)^ 


vWf; 


-(WfGbD)  faOi’V)^) +swf2(V,,p)(^)j]l 
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COMPUTER  PROGRAM  'FINIMP' 
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